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^ ■ Abstract 

o 

' In very high energy nuclear cohisions, the initial energy of produced gluons per 

O . unit area per unit rapidity, dE/L'^/drj, is equal to f {g"^ iJ,L) {g"^ / g"^ , where //^ 

is proportional to the gluon density per unit area of the colliding nuclei. For an 
SU(2) gauge theory, we perform a non-perturbative numerical computation of the 
p • function /{g'^fiL). It decreases rapidly for small g^fiL but varies only by ~ 25%, 

jL^: from 0.208±0.004 to 0.257±0.005, for a wide range 35.36-296.98 in g'^^L, including 

I the range relevant for collisions at RHIC and LHC. Extrapolating to SU(3), we 

>-Ch ■ estimate the initial energy per unit rapidity for Au-Au collisions in the central 

region at RHIC and LHC. 

X 

H 

\ By the end of 1999, the Relativistic Heavy Ion collider (RHIC) at BNL will begin 

colliding beams of Au-ions at ^/s = 200 GeV/nucleon. Some years later, the Large 
Hadron collider (LHC) at CERN will collide heavy ions at ^/s ~ 5.5 TeV/nucleon. The 
objective of these experiments is to understand the properties of very hot and dense 
partonic matter in QCD. It is of considerable interest to determine whether this hot and 
dense matter equilibrates to briefly form a plasma of quarks and gluons (QGP) 

The dynamical evolution of such a system clearly depends on the initial conditions, 
namely, the parton distributions in the nuclei prior to the collision. For partons with 
transverse momenta pt ^ ^qcd, cross-sections in the standard perturbative QCD ap- 
proach may be computed by convolving the parton distributions of the two nuclei with 
the elementary parton-parton scattering cross sections. At the high energies of the RHIC 
(LHC) collider, hundreds (thousands) of mini-jets with pts of the order of several GeV 
may be formed [0]. Final state interactions of these mini-jets are often described in multi- 
ple scattering (see Ref. and references therein) or in classical cascade approaches (see 
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Ref. and references therein). Estimates for the initial energy density in a self screened 
parton cascade approach can be found in Ref. 0. 

At central rapidities, where x <^ 1, and pt ^ ^qcd with x defined to be Pt/i/i, 
parton distributions grow rapidly, and may even saturate for large nuclei for x's in the 
range 10~^ to 10~^ relevant for nuclear collisions at RHIC and LHC respectively 0. 
Coherence effects are important here, and are only included heuristically in the above 
mentioned perturbative approaches. 

In this letter, we will describe results from a classical effective field theory (EFT) 
approach which includes coherent effects in the small x parton distributions of large 
nuclei If the parton density in the colliding nuclei is large at small x, classical meth- 
ods are applicable. It has been shown recently that a renormalization group improved 
generalization of this effective action reproduces several key results in small x QCD: the 
leading a5'log(l/x) BFKL equation, the double log GLR equation and its extensions, 
and the small x DGLAP equation for quark distributions P, |lOl- It has also been argued. 



from other considerations, that the main results of this model should be general results 
in small x QCD g]. 

The EFT contains one dimensionful parameter /z^, which is the variance of a Gaussian 
weight over the color charges of partons, of each nucleus, at rapidities higher than 



the rapidity of interest. For central impact parameters, it is determined to be |^ 




q{x,Q') + j^^gM')^ , (1) 

where Q"^) and xg{x, Q'^) stand for the nucleon quark and gluon structure functions 
at the resolution scale Q of the physical process of interest. Also, one has xq = Q/y/s, 
tq = 1.12 fm, and Nc is the number of colors. From the HERA data for g & one 
obtains /i < 1 GeV for LHC energies and /i < 0.5 GeV at RHIC [|11|. The classical gauge 
fields, and hence the classical parton distributions, can be determined analytically |12| . 
On this basis, it has been argued recently that the typical transverse momenta scale Qg 
in this model is further in the weak coupling regime, with Qg ^ 1 GeV for RHIC and 
Qs ~ 2-3 GeV at LHC 



Kovner, McLerran and Weigert [|T^ applied the effective action approach to nuclear 



collisions. (For an interesting alternative approach, see Ref. [T^]-) Assuming boost 



invariance, and matching the equations of motion in the forward and backward light 
cone, they obtained the following initial conditions for the gauge fields in the = 
gauge: ^4^1^=0 = A\ + Ai , and v4±|^=o = ± f a;±[v4i, A^]. Here A*i,2(p^) {i = 1, 2) are the 
pure gauge transverse gauge fields corresponding to small x modes of incoming nuclei 
(with light cone sources Slx"^)) in the 6{±x^)6{^x^) regions respectively of the light 
cone. 

The sum of two pure gauges in QCD is not a pure gauge-the initial conditions there- 
fore give rise to classical gluon radiation in the forward light cone. For pt » asfJ^, 
the Yang-Mills equations may be solved perturbatively to quadratic order in asfJ^/pt- 
After averaging over the Gaussian random sources of color charge on the light cone. 
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the perturbative energy and number distributions of physical gluons were computed by 



several authors |13|, |TT], |T5|, |16|, 0. In the small x limit, it was shown that the clas- 
sical Yang-Mills result agreed with the quantum Bremsstrahlung result of Gunion and 
Bertsch [||]. 

In Ref. [|l^, we suggested a lattice discretization of the classical EFT, suitable for 
a non-perturbative numerical solution. Assuming boost invariance, we showed that in 
A'^ = gauge, the real time evolution of the small x gauge fields A^{xt,T), A^^XtyT) is 
described by the Kogut-Susskind Hamiltonian in 2-|-l-dimensions coupled to an adjoint 
scalar field. The lattice equations of motion for the fields are then determined straight- 
forwardly by computing the Poisson brackets. The initial conditions for the evolution 
are provided by the lattice analogue of the continuum relations discussed earlier in the 
text. We impose periodic boundary conditions on an x transverse lattice, where A^ 
denotes the number of sites. The physical linear size of the system is L = a A^, where a 
is the lattice spacing. It was shown in Ref. that numerical computations on a trans- 



verse lattice agreed with lattice perturbation theory at large transverse momentum. For 
details of the numerical procedure, and other details, we refer the reader to Ref. [pO[] . 

In this letter, we will focus on computing the energy density e as a function of the 
proper time r. This computation on the lattice is straightforward. Our main result 
is contained in Eq. To obtain this result, we compute the Hamiltonian density on 
the lattice for each p^, and then take the Gaussian average (with the weight /i^) over 
between 40 p trajectories for the larger lattices and 160 p trajectories for the smallest 
ones. 

In our numerical simulations, all the relevant physical information is compressed in 
g^p, and L, and in their dimensionless product g^pL ||2l|. The strong coupling constant 



g depends on the hard scale of interest; from Eq. |I], we see that p depends on the nuclear 
size, the center of mass energy, and the hard scale of interest; is the transverse area 
of the nucleus fsg. Assuming ^ = 2 (or = l/vr), p = 0.5 GeV (1.0 GeV) for RHIC 



(LHC), and L = 11.6 fm for Au-nudei, we find g'^pL ^ 120 for RHIC and ^ 240 
for LHC. (The latter number would be smaller for a smaller value of g at the typical 
LHC momentum scale.) As will be discussed later, these values of g^pL correspond to 
a region in which one expects large non-perturbative contributions from a sum to all 
orders in ~ 6asp/pt, even if 0:5 <^ 1. We should mention here that deviations from 
lattice perturbation theory, as a function of increasing g^pL, were observed in our earlier 
work EOi. 



We shall now discuss results from our numerical simulations. In Fig. 1, we plot 
er/{g^p)^, as a function of g'^pr, in dimensionless units, for the smallest, largest, and 
an intermediate value in the range of g^pUs studied. The quantity er has the physical 
interpretation of the energy density of produced gluons dE/L?'/dri only at late times- 
when T t. Though er goes to a constant in all three cases, the approach to the 
asymptotic value is different. For the smallest g^pL, er increases continously before 
saturating at late times. For larger values of g^pL, er increases rapidly, develops a 
transient peak at r ~ l/g^p, and decays exponentially there onwards, satisfying the 
relation a + /^e""^^, to a constant value a (equal to the lattice dE / L'^ / drj\) . The lines 
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shown in the figure are from an exponential fit including all the points past the peak. 
This behavior is satisfied for all g^^L > 8.84, independently of A^. 

Given the excellent exponential fit, one can interpret the decay time r^) = l/j/g'^fi as 
the appropriate scale controlling the formation of gluons with a physically well defined 
energy. In other words, is the "formation time" in the sense used by Bjorken |2^, ^| 



In Table 1, we tabulate 7 versus g^fiL for the largest N x N lattices for all but the 
smallest g^^L. For large g^^iL, the formation time decreases with increasing g'^fiL, as 
we expect it should. 

In Fig. 2, we plot the asymptotic values a oier / {g'^[j)^ as a function of g'^^ia for various 
values of g^nL. As shown in the upper part of Fig. 2, for smaller g'^fiL, one can go very 
close to the continuum limit with excellent statistics (over 160 independent p trajectories 
for the two smallest values of g'^fiL). In the lower part of Fig. 2, all the data give straight 
line fits with good x^squared values. We use these fits to extrapolate the value of a in 
the continuum limit. We note that the largest value of g^pL with the smallest g^pa 
equal to 0.247, is relatively much further away from the continuum limit than the points 
in the upper part of the figure. It is obtained by averaging 40 independent trajectories 
on a 1200 x 1200 lattice. To lower g'^fia below 0.1, would require going to lattices with 
3000 X 3000 sites. This exceeds the CPU memory of our current computational resources. 
Nevertheless, even for the largest g'^pL, we do get a fine linear fit-though we would warn 
of a potentially large systematic error in the extrapolated value of 6T/{g'^fi)^. 

The physical energy per unit area per unit rapidity of produced gluons can be defined 
in terms of a function f{g'^fiL) as 

The function / here is obtained by extrapolating the values in Fig. 2 to the continuum 
limit. In Fig. 3, we plot the striking behavior of / with g^pL. For very small g'^fiL's, 
it changes very slightly but then changes rapidly by a factor of two from 0.427 to 0.208 
when g^pL is changed from 8.84 to 35.36. From 35.36 to 296.98, nearly an order of 
magnitude in g'^fiL, it changes by ~ 25%. The precise values of / and the errors are 
tabulated in Table 1. 



g'^fiL 


5.66 




8.84 


17.68 


35.36 




70.7 


f 


.436 ± .007 


.427 ± 


004 


.323 ± .004 


.208 ± .004 


.200 ± 


.005 


7 




.101 ± 


024 


.232 ± .046 


.165 ± .013 


.275 ± 


.011 


g'pL 


106.06 


14 


8.49 


212.13 


296.98 






f 


.211 ± .001 


.232 ± 


001 


.234 ± .002 


.257 ± .005 






7 


.322 ± .012 


.362 ± 


023 


.375 ± .038 


.378 ± .053 







Table 1: The function / = dE/L'^/df] and the relaxation rate 7 = l/Tjy/g'^n tabulated as 
a function of g^pL. 7 has no entry for the smallest g^fiL since there eT/{g'^pY g^pr 
differs qualitatively from the other g^^L values. 
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What is responsible for the dramatic change in the behavior of / as a function of 
g^liLl In A"^ = gauge, the dynamical evolution of the gauge fields depends entirely 
on the initial conditions, namely, the parton distributions in the wavefunctions of the 
incoming nuclei [Q. In the nuclear wavefunction, at small x, non-perturbative, albeit 



weak coupling, effects become important for transverse momenta ~ 6 a^/i. The EFT 
predicts that classical parton distributions which have the characteristic Weizsacker- 
Williams l/p^ behaviour for large transverse momenta [pt 3> Qs) grow only logarithmi- 
cally for pt < Qs- One can therefore think of Qs as a saturation scale p9[ that tempers 
the growth of parton distributions at small momenta. 

Now on the lattice, pt is defined to be 2iTn/L, where n labels the momentum mode. 
The condition that momenta in the wavefunctions of the incoming nuclei have saturated, 
Pt ~ 6 05/2, translates roughly into the requirement that g^^L > 13 for n = 1. Thus for 
g^nL = 13, one is only beginning to sample those modes. Indeed, this is the region in 
g^nL in which one sees the rapid change in /. The rapid decrease in / is likely because 
the first non-perturbative corrections are large, and have a negative sign relative to 
the leading term. Understanding the later slow rise and apparent saturation with g^fiL 
requires a better understanding of the number and energy distributions with pt. This 
work is in progress and will be reported on separately 



Our results are consistent with an estimate by A. H. Mueller |^5| for the number 
of produced gluons per unit area per unit rapidity. He obtains dN/L'^/dr] = c{N^ — 
1) Q'^/An'^ as Nc, and argues that the number c is a non-perturbative constant of order 
unity. If most of the gluons have pt ~ Qs, then dE/L'^/d-q = c' {N^ — 1) Q^/Atc^ as Nc 
which is of the same form as our Eq. 0. In the g^^L region of interest, our function 
/ ^ 0.23-0.26. Using the relation between Qs and g^fi [29|, we obtain c' = 4.3-4.9. 
Since one expects a distribution in momenta about Qs, it is very likely that c' is at least 
a factor of 2 greater than c-thereby yielding a number of order unity for c as estimated 
by Mueller. This coefficient can be determined more precisely when we compute the 
non-perturbative number and energy distributions. 

We will now estimate the initial energy per unit rapidity of produced gluons at 
RHIC and LHC energies. We do so by extrapolating from our SU(2) results to SU(3) 
assuming the Nc dependence to be {N^ — 1)/Nc as in Mueller's formula. At late times, 
the energy density is £ = (g'^fj,)^ f{g^fiL) ■j{g'^fiL)/g'^, where the formation time is td = 
l/^{g'^liL) / g'^li as discussed earlier. We find that e^^^^ ^ 66.49 GeV/fm^ and e^^^ « 
1315.56 GeV/fm^. Multiplying these numbers by the initial volumes at the formation 
time Td, we obtain the classical Yang-Mills estimate for the initial energies per unit 
rapidity Et to be E^^^^ ^ 2703 GeV and E!^"^ ^ 24572 GeV respectively. 



Compare these numbers to results presented recently by Kajantie |]26| for the mini- 
jet energy (computed for pt > Psat, where Psat is a saturation scale akin to Qs). He 
obtains E^^^'-^ = 2500 GeV and E^f^'-' = 12000. The remarkable closeness between our 
results for RHIC is very likely a coincidence. Kajantie's result includes a K factor of 
1.5-estimates range from 1.5-2.5 [^. If we pick a recent value of ~ 2 [^, we obtain 
as our final estimate, E^^^^ ^ 5406 GeV and E^^^ ^ 49144 GeV. 
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To summarize, we performed a non-perturbative, numerical computation, for a SU(2) 
gauge theory, of the initial energy, per unit rapidity, of gluons produced in very high 
energy nuclear collisions. Extrapolating our results to SU(3), we estimated the ini- 
tial energy per unit rapidity at RHIC and LHC We plan to improve our estimates by 
performing our numerical analysis for SU(3). Moreover, computations in progress to 
determine the energy and number distributions should enable us to match our results at 
large transverse momenta to mini-jet calculations P^|. 
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Figure 1: er/^g'^fi)^ as a function of g'^fir for g'^fiL = 5.66 (diamonds), 35.36 (pluses) 
and 296.98 (squares). Both axes are in dimensionless units. Note that £r = at r = 
for all g'^/jiL. The lines are exponential fits a + P e"^'^ including all points beyond the 
peak. 
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Figure 2: eT/{g'^fi)^ as a function of g'^fia. The points in the upper plot correspond to 
g^fiL = 5.66 (diamonds), 8.84 (pluses), 17.68 (squares), and 35.36 (x's). The lower plot 
has g^fiL = 70.7 (diamonds), 106.06 (pluses), 148.49 (squares), 212.13 (triangles) and 
296.98 (x's). Lines in the lower plot are fits of form a — b-x. The g'^fJ-a ranges are different 
in the two halves. The points in the upper half are typically closer to the continuum 
limit. 
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Figure 3: erj^g^iif' extrapolated to the continuum limit: / as a function of g^^iL. The 
error bars are smaller than the plotting symbols. 
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